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Hysteresis  effects  and  reaction  characteristics  of  methane  catalytic  partial  oxidation  (CPO)  in  a  fixed-bed 
reactor  are  numerically  simulated.  The  reactions  are  modeled  based  on  the  experimental  measurements 
of  methane  CPO  with  a  rhodium  (Rh)  catalyst.  Three  C/O  ratios  of  0.6,  1.0  and  1.4  are  considered  in 
the  study.  When  the  Reynolds  number  is  200,  the  predictions  indicate  that  the  methane  CPO  is  always 
triggered  at  around  the  inlet  temperature  of  550  K,  regardless  of  what  the  C/O  ratio  is.  It  is  of  interest 
that  if  the  inlet  temperature  is  decreased  after  the  methane  CPO  develops  at  higher  inlet  temperatures, 
the  reversed  path  of  methane  conversion  is  different  from  the  original  path  at  lower  inlet  temperatures. 
The  hysteresis  effect  of  the  methane  CPO  is  thus  observed.  The  hysteresis  behavior  implies  that  a  higher 
yield  of  syngas  or  hydrogen  can  be  achieved  by  controlling  the  reaction  process.  Decreasing  the  C/O  ratio 
intensifies  the  methane  CPO  so  that  the  hysteresis  effect  is  more  pronounced,  and  vice  versa.  An  increase 
in  Reynolds  number  delays  the  excitation  temperature  of  methane  CPO  and  lessens  the  hysteresis  effect 
of  methane  conversion  due  to  the  growth  of  fluid  inertial  force.  However,  the  hysteresis  effect  of  the 
maximum  temperature  in  the  catalyst  bed  increases  as  a  result  of  more  methane  consumption. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Methane  is  a  common  fuel  which  has  been  widely  used  for 
the  purpose  of  getting  heat,  power  and  electricity.  Methane  is  also 
an  important  raw  material  in  the  chemical  industry.  For  example, 
from  the  steam  reforming  of  methane,  synthesis  gas  or  syngas, 
a  mixture  of  hydrogen  and  carbon  monoxide,  can  be  produced 
[1].  Syngas  is  commonly  used  for  methanol  and  ethanol  syn¬ 
thesis  as  well  as  Fisher-Tropsch  process  [2,3].  In  recent  years, 
because  of  the  increasing  interest  in  a  more  efficient  and  low  pol¬ 
lution  energy  conversion  system,  there  has  been  a  pronounced 
progress  in  the  development  of  fuel  cells.  Fuel  cells  can  be  oper¬ 
ated  in  either  low-,  medium-  or  high-temperature  environments, 
depending  on  the  types  of  fuel  cell  [4].  For  high-temperature 
fuel  cells  such  as  molten  carbonate  fuel  cells  (MCFCs)  and  solid 
oxide  fuel  cells  (SOFCs),  syngas  can  act  as  a  feedstock  to  drive 
the  devices  and  generate  electricity  [5].  Syngas  can  further  be 
processed  through  water  gas  shift  reactions  [6,7]  to  promote  the 
concentration  of  hydrogen.  Once  the  hydrogen  has  been  separated 
and  purified,  it  can  be  used  as  a  fuel  in  low-  and  medium- 
temperature  fuel  cells  such  as  a  proton  exchange  membrane  fuel  cell 
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(PEMFC),  an  alkaline  fuel  cell  (AFC)  and  a  phosphoric  acid  fuel  cell 
(PAFC). 

Though  steam  methane  reforming  (SMR)  is  an  important  tech¬ 
nique  for  producing  syngas  and  hydrogen,  the  reaction  pertains  to 
a  highly  endothermic  reaction  in  nature  [8]. 

CH4  +  H2O  CO  +  3H2,  AH  =  +206  kj mol”1  (1) 

On  account  of  the  strong  endothermic  reaction  of  SMR,  extra 
heat  has  to  be  supplied  while  the  syngas  is  generated.  To  overcome 
this  drawback,  syngas  or  hydrogen  generation  from  methane  can 
be  carried  out  by  means  of  the  method  of  catalytic  partial  oxidation 
(CPO)  [9,10].  The  process  of  methane  CPO  is 

CH4  +  Io2  -*  CO  +  2H2,  AH  =  36kJmol_1  (2) 

As  a  whole,  the  process  is  a  mild  exothermic  reaction.  In  other 
words,  instead  of  heat  being  supplied  from  a  separate  energy 
source,  methane  CPO  is  a  naturally  autothermal  process.  When  the 
methane  CPO  is  further  analyzed,  the  detailed  reactions  include 
methane  combustion,  steam  reforming  and  carbon  dioxide  reform¬ 
ing  (or  dry  reforming)  [8].  Specifically,  the  energy  required  for 
producing  syngas  is  provided  from  the  partial  oxidation  or  com¬ 
bustion  of  methane.  Methane  combustion  is  a  highly  exothermic 
reaction  and  it  is  expressed  as 

CH4+202  ->•  C02  +2H20,  AH  =  -803  kj  mol-1  (3) 
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Nomenclature 

Dj  diffusion  coefficient  of  species  i  (m2  s-1 ) 

Dp  mean  particle  diameter  (m) 

E  total  energy  (J) 

Fa  activation  energy  (kj  mol-1 ) 

Fj  mole  fraction 

hi  enthalpy  of  species  i  (J  kg-1 ) 

standard-state  enthalpy  of  species  i  (J  kg-1 ) 

AH0  change  of  total  enthalpy  on  standard  state 
Ji  diffusion  flux  of  species  i  (kg  m-2  s-1 ) 

I<e ff  effective  thermal  conductivity  (W  m-1  I<-1 ) 
kf  fluid  phase  thermal  conductivity  (W  m-1 I<-1 ) 

ks  solid  medium  thermal  conductivity  (W  m-1  K-1 ) 

/<eq  equilibrium  constant 

Mi  molecular  weight  of  species  i  (kg  kgmol-1 ) 

N  number  of  species 

p  pressure  (Pa) 

r  radial  coordinate  (m) 

R  universal  gas  constant  (8.314  m3  Pa  K-1  mol-1 ) 

Ri  reaction  rate  of  species  i  (kmol  m-3  s-1 ) 

S  selectivity 

AS0  change  of  total  entropy  on  standard  state 

T  temperature  (K) 

v  velocity  (ms-1) 

x  axial  coordinate  (m) 

Y;  mass  fraction  of  species  i 

Greek  letters 
y  porosity 

p  density  (kg  m-3) 

r  shear  stress  (Pa) 

fi  viscosity  (Pas) 

Subscript 

CH4  methane 

CO  carbon  monoxide 

com  combustion 

cr  carbon  dioxide  reforming 

f  fluid 

H20  water 

i  species  i 

in  inlet 

N2  nitrogen 

out  outlet 

sr  steam  reforming 

steam  steam 

w  wall 


With  the  aid  of  catalysts,  syngas  formation  from  methane  can  be 
triggered  at  a  lower  temperature  environment  and  the  heat  released 
from  the  combustion  is  able  to  induce  the  subsequently  endother¬ 
mic  reactions  of  steam  reforming  (i.e.  Eq.  (1))  and  carbon  dioxide 
reforming.  The  reaction  of  carbon  dioxide  reforming  is  written  as 
the  following 

CH4  +  C02^>  2 CO  +  2H2,  AH  =  247  kj  mol-1  (4) 

Accordingly,  it  is  known  that  the  products  of  methane  combus¬ 
tion  (i.e.  H20  and  C02)  are  then  reduced  to  hydrogen  and  carbon 
monoxide,  thereby  yielding  syngas. 

It  is  known  that  catalysts,  which  are  used  in  methane  CPO,  play 
a  key  role  in  determining  the  methane  conversion  and  syngas  for¬ 
mation.  In  this  aspect,  some  group  VIII  or  noble  metals  such  as  Ni, 
Pt,  Ru  and  Rh  have  been  widely  employed  and  they  are  usually  sup¬ 


ported  on  certain  oxides  such  as  a-Al203  and  Ru02.  For  instance, 
in  the  earlier  of  study  of  Prettre  et  al.  [11  ],  methane  CPO  was  con¬ 
ducted  over  10  wt%  refractory  supported  nickel  in  the  temperature 
range  of  973-1173  K  at  1  atm.  Ashcroft  et  al.  [12]  studied  methane 
CPO  with  the  catalysts  of  lanthanide  ruthenium  oxides  as  well  as 
Ru  supported  on  A1203  and  pure  Ru02  at  1050  K  and  1  atm.  They 
reported  that  methane  conversion  was  in  excess  of  90%  and  syngas 
selectivity  was  in  the  range  of  94-99%.  In  the  experimental  study 
of  Horn  et  al.  [13],  methane  CPO  on  Rh  and  Pt  foam  catalysts  were 
analyzed  and  compared  with  each  other.  They  reported  that  Rh  cat¬ 
alyzed  steam  reforming  much  more  efficiently  than  Pt  so  that  the 
hydrogen  and  carbon  monoxide  yields  were  generally  higher  on  the 
former  than  on  the  latter. 

In  addition  to  experimental  studies,  a  number  of  numerical 
simulations  have  also  been  performed  to  figure  out  the  detailed 
reaction  processes  of  methane  CPO.  In  the  earlier  studies  of  Hick¬ 
man  and  Schmidt  [  14,1 5  ] ,  methane  CPO  over  Pt  and  Rh  surfaces  were 
predicted  using  a  model  which  consisted  of  19  elementary  steps 
including  adsorption,  desorption  and  surface  reactions.  It  was  illus¬ 
trated  that,  under  methane-rich  conditions  at  high-temperatures, 
H2  and  CO  were  the  primary  products  of  methane  direct  oxidation 
through  a  pyrolysis  mechanism.  Potential  energy  diagrams  were 
also  displayed  in  their  studies.  The  study  of  Williams  et  al.  [16] 
further  developed  two  different  mechanisms  to  predict  methane 
CPO  on  Rh-coated  foam  monoliths;  one  was  a  38-step  mechanism 
and  the  other  a  104-step  mechanism.  They  highlighted  the  impor¬ 
tance  of  high-resolution  spatial  profiles  for  each  mechanism.  Unlike 
the  aforementioned  multi-step  mechanisms  developed,  Jin  et  al. 
[17]  modeled  and  simulated  methane  CPO  in  a  membrane  reactor 
using  one-step  global  chemical  mechanisms.  In  spite  of  such  sim¬ 
ple  kinetics,  their  simulated  results  were  in  good  agreement  with 
the  experimental  data  and  it  was  found  that  at  a  low  carbon  space 
velocity  and  when  the  temperature  was  between  1098  and  1158  K, 
CH4  conversion  was  larger  than  96%  and  CO  selectivity  was  larger 
than  97%.  Abashar  et  al.  [18]  used  the  developed  model  of  Jin  et  al. 
[17]  to  predict  the  coupled  steam  and  oxidative  reforming  reactions 
of  methane  in  a  circulating,  fast  fluidized  bed  membrane  reactor 
where  an  oxygen  discrete  injection  was  proposed  to  overcome  the 
hot  spot  temperature  problem. 

The  above  literature  has  provided  some  important  information 
to  describe  the  reaction  characteristics  of  methane  CPO  and  these 
studies  are  conducive  to  practicing  the  production  of  syngas  and 
hydrogen  for  industrial  purposes.  As  illustrated  above,  Rh  catalysts 
are  usually  more  appropriate  than  Pt  catalysts  for  syngas  gener¬ 
ation  [13].  As  a  result,  the  present  study  is  intended  to  establish 
reaction  models  in  terms  of  a  number  of  one-step  global  chemical 
mechanisms  to  describe  methane  CPO  on  Rh.  Particular  emphasis 
is  placed  on  the  hysteresis  effects  of  methane  CPO  at  various  inlet 
conditions.  Details  of  the  hysteresis  effect  and  reaction  behavior  of 
the  methane  CPO  on  Rh  will  be  pointed  out. 

2.  Mathematical  formulation  and  modeling 

2.1.  Geometry  and  assumptions 

A  schematic  of  the  computational  geometry  and  physical  size 
of  the  studied  fixed-bed  reactor  is  sketched  in  Fig.  la  where  the 
simulation  is  performed  in  a  cylindrical  tube.  The  geometry  of  the 
reactor  is  constructed  based  on  the  experimental  set-up  of  Horn  et 
al.  [13].  As  a  whole,  the  reactor  consists  of  three  zones.  The  middle 
(second)  zone  is  filled  with  the  Rh  catalyst  which  is  sandwiched  by 
the  first  and  third  zones.  The  first  and  third  zones  remain  empty 
for  fluid  to  flow  inward  and  outward  from  the  catalyst  bed,  respec¬ 
tively.  These  three  zones  are  enveloped  by  a  cylindrical  shell  with 
an  adiabatic  wall.  Therefore,  the  released  heat  will  not  penetrate 
through  the  tube  wall. 
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Fig.  1.  Schematics  of  (a)  computational  geometry  and  reactor  size  as  well  as  (b) 
adopted  grid  system  (120  x  33). 


Seeing  that  the  investigated  phenomena  are  related  to  fluid 
motion,  heat  and  mass  transfers  as  well  as  chemical  reactions,  some 
assumptions  are  adopted  to  simplify  the  physical  problem.  They 
include:  (1)  the  flow  field  is  laminar  and  axisymmetric  along  the 
centerline  of  the  reactor,  whereby  the  two-dimensional  simulation 
is  carried  out;  (2)  the  porous  media  in  the  catalyst  bed  is  homoge¬ 
neous  and  thermal  equilibrium  prevails  at  the  catalyst  surface  so 
that  there  is  no  temperature  difference  at  the  interface  of  the  cata¬ 
lyst  and  the  fluid  [4,13  ] ;  (3 )  the  effect  of  gravity  on  the  flow  field  and 
the  reactions  of  reagents  are  ignored;  and  (4)  gas  mixture  inside  the 
reactor  abides  by  the  ideal  gas  law. 

2.2.  Governing  equations 


Energy  equation: 

N  h° 

V  •  (v(pfEf  +  p))  =  V  •  [keffVT  -  YhJi)  +  ■ t)]  -  J2wRi  (10) 

i  i 

where  I<eff=ykf+{\ -y)I<s.  In  both  the  non-porous  and  porous 
regions,  the  species  equation  and  the  equation  of  state  are  also 
invoked.  They  are  given  as  follows. 

Species  equation: 

N 

V  •  (jivYi)  =  V  ■  (pDjVYj)  +  M^Ri  (U ) 

Z  =  1 

Equation  of  state: 

(,2> 

i 

2.3.  Boundary  conditions 


When  boundary  conditions  are  considered,  four  different  sec¬ 
tions  are  involved,  including  the  upstream  inflow,  the  downstream 
outflow,  the  centerline  and  the  wall.  For  the  upstream  inflow,  the 
boundary  conditions  are  given  by 


u  =  uin,v  =  vjn,Yi  =  YiM,T  =  Tin,  and  p  =  pin  (13) 

In  regard  to  the  downstream  outflow,  they  are 
V  -v  =  VYj  =  VT  =  0  (14) 


By  virtue  of  the  absence  of  heat  and  mass  fluxes  across  the  cen¬ 
terline,  the  boundary  conditions  along  the  axissymmetric  axis  can 
be  treated  as  the  following 


dT 

dr 


(15) 


Inside  the  reactor  surface,  the  no-slip  velocity  and  adiabatic  wall 
are  encountered.  Hence,  the  boundary  conditions  are  expressed  to 
be 


dYi  dT 


(16) 


By  virtue  of  the  chemical  reactions  inside  a  catalyst  bed  encoun¬ 
tered,  two  kinds  of  governing  equations  are  considered  to  take 
the  non-porous  and  porous  zones  into  account.  In  the  non-porous 
regions,  the  governing  equations  include 
Continuity  equation: 

V.(pv)  =  0  (5) 

Momentum  equation: 

V  •  (pvv)  =  -Vp  +  V  •  (t)  (6) 

Energy  equation: 

V  .  (v(pE  +p))  =  V  .  [kVT  -  (5>)  +  (f  •  i))  (7) 


In  the  porous  region  or  the  catalyst  bed,  the  governing  equations 
contain 

Continuity  equation: 

V-(yp5)  =  0  (8) 


Momentum  equation: 


150  •  p  (1  -  y)2 
Dj  Y3 


Poo 


(9) 


2.4.  Model  of  methane  CPO 

As  far  as  the  chemical  reactions  are  concerned,  it  is  known  that 
accurate  predictions  of  physical  phenomena  are  highly  related  to 
reaction  rate  expressions  and  the  chemical  reaction  rates  depend 
strongly  on  the  catalyst  employed.  Presently,  numerical  simulations 
of  methane  CPO  on  an  Rh  catalyst  are  performed.  Similar  to  the 
one-step  global  chemical  mechanisms  adopted  by  Jin  et  al.  [17]  and 
Abashar  et  al.  [  18  ] ,  the  reaction  rates  of  methane  combustion  (RCom ), 
steam  reforming  (RSr)  and  carbon  dioxide  reforming  (RCr)  are  also 
modeled  using  one-step  global  chemical  mechanisms.  However,  the 
rate  expressions  are  established  based  on  the  experimental  data  of 
Horn  et  al.  [13]  and  they  are  given  in  the  following  forms 

Rcom  =  2.119  X  1010exp  (2^)  Pch4Po2  (17) 

Rsr  =  1.0  x  10“9r5exp  (-^)  (PCH4PH20  -  Ke~q!srPCH4PH2o) 

(18) 

Rcr  =  7.0  x  10-7r4exp  (--^)  (PCH4PC02  -  ^’crPCH4PC02) 


V  ■  ( ypvv )  =  -yVp  +  V  •  (y?) 


(19) 
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Fig.  2.  Distributions  of  methane  conversion  along  the  centerline  of  the  catalyst  bed 
at  three  different  grid  systems. 

where  JCeq  is  the  equilibrium  constant  which  links  the  forward 
and  backward  reactions.  The  constant  is  obtained  in  accordance 
with  the  thermodynamic  theory  [19]  and  it  is  expressed  as 

2.5.  Numerical  method 

To  predict  the  methane  CPO  phenomena  in  association  with  the 
aforementioned  models,  the  commercial  software  FLUENT  (Ver¬ 
sion  6.2)  is  utilized  to  solve  the  governing  equations  along  with  the 
boundary  conditions.  In  brief,  the  finite-volume  method  in  associa¬ 
tion  with  SIMPLE  algorithm  [20]  is  utilized.  The  first  order  upwind 
scheme  is  used  to  evaluate  the  convective  and  diffusive  fluxes  over 
the  control  volume  surface.  An  orthogonal  grid  system  is  employed 
to  reduce  numerical  truncation  errors,  which  may  occur  during 
the  course  of  calculating.  Regarding  the  grid  system,  the  grids  are 
uniformly  distributed  along  the  x-  and  y-axes.  Three  different  grid 
systems  of  60  x  17, 120  x  33  and  240  x  66  are  tested  and  compared 
with  each  other  to  seek  a  proper  grid  system  for  simulations.  With 
the  three  grid  systems,  the  distributions  of  methane  conversion 
along  the  centerline  of  the  reactor  in  the  catalyst  bed  are  shown 
in  Fig.  2  where  the  inlet  temperature  is  300  °C  and  the  C/O  ratio  is  1. 
It  is  obvious  that  no  matter  which  grid  system  is  adopted,  the  dis¬ 
crepancy  among  the  three  curves  is  almost  imperceptible.  In  view  of 
the  grid  system  of  120  x  33  satisfying  the  requirement  of  grid  inde¬ 
pendence,  it  is  employed  to  simulate  the  phenomena  of  methane 
CPO. 


Fig.  3.  Comparisons  of  (a)  methane  conversion  and  (b)  CO  selectivity  between  the 
numerical  simulation  and  the  experimental  measurement. 

To  validate  the  developed  models  and  adopted  numerical  meth¬ 
ods,  comparisons  between  the  numerical  predictions  and  the 
experimental  measurements  [13]  are  displayed  in  Fig.  3  where  the 
methane  conversion  and  CO  selectivity  at  the  exit  of  the  reactor 
are  plotted.  In  the  experiments  of  Florn  et  al.  [13],  methane  mixed 
with  oxygen  and  argon  passing  through  the  Rh  catalyst  bed  was 
carried  out.  As  a  whole,  the  numerical  predictions  of  the  methane 
conversion  are  in  agreement  with  the  experimental  data  (Fig.  3a). 


Table  1 

Inlet  conditions  of  methane  CPO  and  the  properties  of  the  catalyst. 


Feed  (based  on  300 1<) 

Feed  flow  rate  (m3  s-1 ) 

Re  GHSV  (h-1 

) 

Feed  composition  (mol%) 

C/O  =  0.6 

C/O  =  1.0 

C/O  =  1.4 

8.91  x  10“5 

200  185,240 

ch4 

20.13 

29.60 

37.03 

02 

16.77 

14.80 

13.22 

n2 

63.10 

55.60 

49.75 

Properties  of  the  catalyst 

Composition  (wt%) 

Porosity  of  catalyst  bed 

Bulk  density  (kg  m-3 ) 

Specific  heat  (J  kg 

r1  K-1) 

Thermal  conductivity  (W  m_  1  K_1 ) 

5%  Rh  +  95%  a-Al203 

0.81 

4325 

848 

36 
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The  predicted  results  of  CO  selectivity  are  somewhat  lower  than 
the  experimental  measurements  (Fig.  3b).  Nevertheless,  the  trend 
of  the  former  is  consistent  with  the  latter.  As  a  consequence,  Fig.  3 
depicts  that  the  established  models  and  adopted  numerical  method 
enable  us  to  investigate  the  phenomena  of  methane  CPO  on  the  Rh 
catalyst  appropriately. 

3.  Results  and  discussion 

Methane  mixed  with  air  flowing  through  the  Rh  catalyst  bed 
serves  as  the  basis  for  the  present  study.  The  stream  in  the  tube  is 
at  atmospheric  pressure  and  the  Reynolds  number  is  fixed  at  200. 
Accordingly,  the  volumetric  flow  rate  of  the  gaseous  mixture  at  the 
inlet  is  8.91  x  10-5  m3  s-1  and  the  gas  hourly  space  velocity  (GFISV) 
is  around  185,240  h-1,  based  on  the  temperature  of  300  K.  Three 
different  C/O  ratios  are  considered;  they  are  0.6, 1.0  and  1.4.  Detailed 
operational  conditions  and  the  properties  of  the  catalyst  are  listed 
in  Table  1. 

3.1  Hysteresis  effect  of  methane  CPO 

Methane  conversion  and  the  maximum  temperature  in  the  reac¬ 
tor  versus  inlet  temperature,  where  the  inlet  temperature  ranges 


from  300  to  800 1<  and  the  C/O  ratio  is  1,  are  sketched  in  Fig.  4a  and 
b,  respectively.  The  methane  conversion  (xCh4)  *s  defined  to  be: 


fo/\  ^CH4,in  ^CH4,out  mn°/ 

XCH4{/°)  =  - F -  x 

FCH4,in 


(21) 


When  the  inlet  temperature  proceeds  from  300  to  500  K, 
the  methane  conversion  is  zero  (Fig.  4a).  This  indicates  that  no 
chemical  reaction  is  triggered,  whereby  the  chemical-frozen  flow 
is  exhibited.  Once  the  inlet  temperature  is  increased  to  550  K, 
methane  CPO  is  triggered  so  that  there  is  a  pronounced  jump  in 


Fig.  4.  Distributions  of  (a)  methane  conversion  and  (b)  the  maximum  temperature  Fig.  5.  Concentration  distributions  of  (a)  H2,  (b)  CO  and  (c)  syngas  versus  inlet 

versus  inlet  temperature  in  the  upper  branch  and  the  lower  branch  at  the  exit.  temperature  in  the  upper  branch  and  the  lower  branch  at  the  exit. 
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methane  conversion.  Specifically,  the  methane  conversion  rises 
from  0  to  69%  (Fig.  4a).  On  the  other  hand,  by  virtue  of  the  par¬ 
tial  combustion  or  violent  oxidation  of  methane  induced,  heat  is 
liberated  from  the  exothermic  reaction  which  results  in  a  sud¬ 
den  growth  in  the  maximum  temperature  (Fig.  4b).  At  the  inlet 
temperature  (550 1<),  the  maximum  temperature  is  around  1590  K. 
Flowever,  this  temperature  is  relatively  low  compared  to  the 
adiabatic  flame  temperature  of  methane.  Typically,  the  constant- 
pressure  adiabatic  flame  temperature  for  the  combustion  of  a 


Inlet  temperature  (K) 


Fig.  6.  Concentration  distributions  of  (a)  CH4,  (b)  CO2  and  (c)  H20  versus  inlet  tem¬ 
perature  in  the  upper  branch  and  the  lower  branch  at  the  exit. 


stoichiometric  methane-air  mixture  is  2318  K  [21  ].  The  lower  max¬ 
imum  temperature  is  attributed  to  the  fact  that  only  part  of  the 
methane  is  depleted  for  combustion,  as  observed  in  Fig.  4a.  Nor¬ 
mally,  the  ignition  temperature  of  methane  is  905  K  [22]  which 
is  much  higher  than  the  activated  temperature  of  the  methane 
CPO  (550  K).  This  indicates  that  the  role,  played  by  the  Rh  cata¬ 
lyst  on  the  methane  CPO,  is  of  the  utmost  importance.  After  the 
methane  CPO  is  excited,  both  the  methane  conversion  and  the 
maximum  temperature  increase  slightly  with  the  increasing  inlet 
temperature. 
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Fig.  7.  Distributions  of  (a)  H2  selectivity,  (b)  CO  selectivity  and  (c)  H2/CO  versus  inlet 
temperature  in  the  upper  branch  and  the  lower  branch  at  the  exit. 
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If  the  inlet  temperature  goes  down  from  800  K,  the  backward 
paths  (the  upper  branch  solutions)  of  both  the  methane  conversion 
and  the  maximum  temperature  are  equivalent  to  the  forward  paths 
(the  lower  branch  solutions).  The  merged  segment  lasts  until  the 
excitation  temperature  of  methane  CPO  (i.e.  550  K).  Once  the  inlet 
temperature  is  lower  than  550  K,  it  is  of  interest  that  the  methane 
conversion  and  the  maximum  temperature  decline  only  a  bit  with 
the  decreasing  inlet  temperature,  even  though  the  inlet  tempera¬ 
ture  goes  as  low  as  300  K.  In  other  words,  the  reversed  paths  are 
totally  different  from  the  original  paths  as  the  inlet  temperature  is 
lower  than  550  K.  Physically,  as  long  as  the  methane  CPO  has  been 
triggered,  the  cooling  effect  brought  by  the  inlet  stream  is  insuf¬ 
ficient  to  eliminate  the  heat  released  from  the  CPO.  As  a  result, 
the  chemical  reaction  sustains  itself.  It  is  thus  recognized  that  the 
methane  CPO  with  respect  to  inlet  temperature  is  characterized  by  a 
hysteresis  effect.  Alternatively,  for  the  presently  studied  conditions 
(i.e.  C/0  =  1  and  Re  =  200),  the  upper  branch  and  the  lower  branch 
are  separated  at  the  excitation  point  of  methane  CPO.  It  follows  that 
the  bifurcation  point  of  the  dual-solution  exactly  stands  for  the  acti¬ 
vated  temperature  of  methane  CPO  in  the  lower  branch.  Reviewing 
the  past  studies  concerning  fuel  burning,  hysteresis  effects  have 
been  observed  in  single-droplet  combustion  [23,24],  two-droplet 
combustion  [25,26],  lifted  flame  [27]  and  swirling  flame  [28]  sys¬ 
tems.  Obviously,  even  though  the  methane  reaction  in  a  CPO  system 
is  not  as  strong  as  those  of  the  aforementioned  burning  systems, 
the  hysteresis  effect  is  still  featured  due  to  the  partial  combustion 
of  methane. 


(b)  Steam  reforming 


X  (m) 


X  (m) 


(c)  C02  reforming 


Reaction  rate 
(kmol  /  m2-s) 

IQ  30 
028 
— 1  026 
□  24 
□  22 
□  20 
□  18 
□  16 
□  14 
□  12 
0.10 
□  □8 
□  □6 
□  □4 

IDD2 
□  □□ 


Reaction  rate 
(kmol  /  m2-s) 


0 


□  24 

□  22 
□  21 
0.19 
□  18 
□  16 

□  14 
0.13 

□  .11 
□  10 
□  □8 
□  □6 
□as 

□  □3 

□  02 
□  oo 


X  (m) 


Fig.  8.  Upper-branching  reaction  rate  contours  of  (a)  methane  combustion,  (b)  steam 
reforming  and  (c)  carbon  dioxide  reforming  in  the  catalyst  bed  with  the  inlet  tem¬ 
perature  of  300  K. 


3.2.  Performance  of  methane  CPO 

Subsequently,  the  exit  concentration  (or  mole  fraction)  distri¬ 
butions  of  H2,  CO  and  syngas  (i.e.  H2  +  CO)  versus  inlet  temperature 
in  the  upper  branch  and  the  lower  branch  are  presented  in  Fig.  5 
to  figure  out  the  performance  of  the  methane  CPO.  When  the  inlet 
temperature  is  lower  than  550  K,  no  reaction  develops  so  that  the 
concentrations  of  H2,  CO  and  syngas  are  zero.  After  the  two  branches 
merge,  increasing  inlet  temperature  makes  these  concentrations 
rise  a  bit.  The  methane  CPO  is  a  slightly  exothermic  reaction,  as 
expressed  in  Eq.  (2).  Basically,  it  is  composed  of  a  forward  com¬ 
bustion  reaction  and  two  reversible  reforming  reactions.  It  should 
be  emphasized  that  the  formation  of  H2  and  CO  comes  from  H20 
and  C02  reforming  (i.e.  Eqs.  (1  and  4)).  The  two  reactions  pertain 
to  an  endothermic  reaction  where  a  higher  reaction  temperature  is 
conducive  to  the  production  of  syngas.  Besides,  the  generation  of 
H20  and  C02  is  subject  to  the  methane  combustion  (i.e.  Eq.  (3)).  In 
other  words,  a  higher  inlet  temperature  facilitates  the  formation  of 
H20  and  C02  in  the  methane  combustion  reaction  and,  furthermore, 
the  steam  and  dry  reforming  reactions  shift  toward  the  formation 
of  H2  and  CO.  This  results  in  the  enhancement  in  the  concentra¬ 
tions  of  H2,  CO  and  syngas  when  the  inlet  temperature  increases. 
The  concentration  distributions  of  CH4,  H20  and  C02  in  the  upper 
branch  and  the  lower  branch  are  provided  in  Fig.  6.  Apparently,  the 
concentration  of  CH4  falls  when  the  methane  CPO  is  excited  in  the 
lower  branch.  On  the  other  hand,  unlike  the  distributions  shown  in 
Fig.  5,  the  concentrations  of  CH4,  H20  and  C02  decline  in  the  upper 
branch  as  the  inlet  temperature  increases,  verifying  that  increasing 
inlet  temperature  enhances  the  forward  reactions  of  H20  and  C02 
reforming. 
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Fig.  9.  Upper-branching  isothermal  contours  in  the  catalyst  bed  with  the  inlet  tem¬ 
peratures  of  (a)  300  K,  (b)  500  K,  and  (c)  700 1<. 
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The  profiles  of  H2  and  CO  selectivity  in  the  two  branches  as  well 
as  H2/CO  ratio  (mole  fraction  ratio)  in  the  upper  branch  are  plotted 
in  Fig.  7.  The  H2  selectivity  (Sh2  )  and  CO  selectivity  (5C0)  are  defined 
as: 

Sh2  =  p — - x  100%  (22) 

^H2,out  +  CH20,out 


Fig.  10.  Concentration  distributions  of  H2,  CO,  C02,  CH4  and  H20  along  the  centerline 
of  the  reactor  in  the  upper  branch  with  the  inlet  temperatures  of  (a)  300  K,  (b)  500  K, 
and  (c)  700  K. 


Sco  =  p - -  x  100%  (23) 

MIO.out  +^C02,out 

It  is  known  that,  following  the  conversion  of  methane,  both 
H2  and  H20  are  generated  when  one  is  concerned  with  hydrogen 
elements.  The  higher  the  H2  selectivity,  more  H20  is  converted 
into  H2.  Within  the  investigated  range  of  the  inlet  temperature 
(300-800  K),  the  H2  selectivity  is  always  higher  than  70%  (Fig.  7a), 
revealing  that  H2  rather  than  H20  is  the  main  product.  It  has  been 
noted  in  Fig.  5  that  a  higher  inlet  temperature  gives  a  higher  max¬ 
imum  temperature,  which  is  conducive  to  the  forward  reaction 
of  H20  reforming  from  the  viewpoint  of  Le  Chatelier’s  principle. 
This  is  the  reason  that  the  H2  selectivity  rises  with  increasing  inlet 
temperature.  The  CO  selectivity  also  grows  with  increasing  inlet 
temperature  and  the  value  is  larger  than  77%  (Fig.  7b),  reflect¬ 
ing  that  CO  formation  dominates  the  competition  between  CO 
and  CO2  generations.  In  examining  Eq.  (2),  it  is  known  that  the 
concentration  ratio  between  FI2  and  CO  from  an  ideal  reaction  of 
methane  partial  oxidation  should  be  2.  Observing  the  curve  shown 
in  Fig.  7c,  an  ascent  in  inlet  temperature  is  closer  to  the  theoretical 
result. 


Fig.  11.  Distributions  of  (a)  methane  conversion  and  (b)  the  maximum  temperature 
versus  inlet  temperature  in  the  upper  and  lower  branch  at  various  C/O  ratios. 
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3.3.  Reaction  characteristic  of  methane  CPO 

To  proceed  farther  into  an  analysis  of  methane  CPO,  the  reaction 
contours  of  methane  combustion,  steam  reforming  and  dry  reform¬ 
ing  in  the  catalyst  bed  with  the  upper  branch  solution  at  300  K  are 
displayed  in  Fig.  8.  Soon  after  the  reactants  flow  into  the  catalyst 
bed,  methane  combustion  occurs  within  a  very  thin  layer,  behaving 
as  a  flame  sheet  (Fig.  8a).  In  contrast,  the  steam  reforming  and  C02 
reforming  prevail  over  the  entire  catalyst  bed.  When  observing  the 
reaction  rate  expressions  in  Eqs.  (17-19),  the  activation  energy  of 
methane  combustion  is  much  higher  than  those  of  steam  and  car¬ 
bon  dioxide  reforming.  It  is  thus  realized  that  all  the  reactions  are 
controlled  by  the  combustion.  Moreover,  on  account  of  inherently 
endothermic  reactions  featured  by  the  steam  reforming  and  dry 
reforming,  only  after  the  combustion  is  triggered,  the  steam  and 
carbon  dioxide  reforming  can  proceed  by  means  of  the  liberated 
heat. 

Fig.  9  shows  the  upper-branching  isothermal  contours  inside  the 
catalyst  bed  at  the  inlet  temperatures  of 300, 500  and  700  K.  Though 
a  lower  inlet  temperature  (300 1<)  causes  a  lower  bed  temperature 
(Fig.  10a)  and  vice  versa  (Fig.  10c),  the  temperature  variation  is  not 


very  significant.  Fig.  10  further  plots  the  concentration  profiles  of 
H2)  CO,  C02,  CH4  and  H20  along  the  centerline  of  the  reactor  at 
the  three  inlet  temperatures.  A  rapid  drop  in  CFI4  concentration  is 
observed  in  the  vicinity  of  the  entrance  of  the  catalyst  bed,  whereas 
the  concentration  distributions  of  FI20  and  C02  are  characterized 
by  a  peak  at  a  certain  location,  as  a  consequence  of  methane  com¬ 
bustion.  Alternatively,  a  progressive  growth  in  the  concentrations 
of  H2  and  CO  is  obtained.  As  a  whole,  the  higher  the  inlet  temper¬ 
ature,  the  more  pronounced  the  concentration  difference  between 
CO  and  C02  as  well  as  H2  and  FI20. 

3.4.  Effects  of  C/O  ratio  and  Reynolds  number  on  hysteresis 
behavior 

The  distributions  of  methane  conversion  and  the  maximum 
temperature  in  the  catalyst  bed  at  various  C/O  ratios  (i.e.  0.6,  1.0 
and  1.4)  are  evaluated  in  Fig.  11,  whereas  the  distributions  of  CO 
and  H2  selectivity  are  displayed  in  Fig.  12.  As  seen  in  Fig.  11  a,  within 
the  investigated  range  of  C/O  ratio,  the  activated  temperature  of 
methane  CPO  is  hardly  affected  by  the  ratio.  Once  the  inlet  tem- 


Fig.  12.  Distributions  of  (a)  CO  selectivity  and  (b)  H2  selectivity  versus  inlet  temper¬ 
ature  in  the  upper  and  lower  branch  at  various  C/O  ratios. 


Fig.  13.  Distributions  of  (a)  methane  conversion  and  (b)  the  maximum  tempera¬ 
ture  versus  inlet  temperature  in  the  upper  and  lower  branch  at  various  Reynolds 
numbers. 
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perature  is  high  enough  to  induce  the  methane  CPO,  the  methane 
conversion  with  C/O  =  0.6  is  much  higher  than  that  with  C/O  =  1.0 
and  1.4,  resulting  from  the  richer  oxygen  contained  in  the  stream.  As 
a  consequence,  the  hysteresis  effect  of  the  methane  CPO  at  C/O  =  0.6 
is  more  profound.  Since  the  methane  CPO  at  C/O  =  0.6  is  closer  to  the 
complete  combustion  of  methane,  it  is  not  surprising  that  the  max¬ 
imum  temperature  of  the  reaction  is  promoted,  in  contrast  to  that 
of  C/O  =  1.0  and  1.4  (Fig.  lib).  Conversely,  from  the  distributions  of 
CO  and  H2  selectivity  shown  in  Fig.  12,  increasing  C/O  ratio  is  con¬ 
ducive  to  the  formation  of  CO  and  H2,  whereby  the  concentration 
of  syngas  in  the  product  gas  rises. 

The  effect  of  varied  Reynolds  number  on  the  hysteresis  phe¬ 
nomena  are  examined  in  Figs.  13  and  14  where  three  different 
Reynolds  numbers  of  200,  800  and  2000  are  considered.  As  the 
Reynolds  number  increases,  the  inertial  force  of  the  inlet  stream 
increases  and  the  residence  time  for  the  reactants  in  the  cata¬ 
lyst  bed  decreases.  As  a  result,  the  methane  CPO  develops  at  a 
higher  inlet  temperature  and  the  hysteresis  behavior  is  less  notable 
when  the  Reynolds  number  is  larger  (Fig.  13a).  Nevertheless,  it  is 
worthy  of  note  that  the  hysteresis  effect  of  the  maximum  temper¬ 
ature  increases  rather  than  decreases  with  increasing  the  Reynolds 


Fig.  14.  Distributions  of  (a)  CO  selectivity  and  (b)  H2  selectivity  versus  inlet  temper¬ 
ature  in  the  upper  and  lower  branch  at  various  Reynolds  numbers. 


number  (Fig.  13b).  This  arises  from  the  fact  that  more  reactants 
are  fed  into  the  catalyst  bed  and  thereby  more  heat  is  liberated 
from  methane  combustion,  even  though  the  methane  conversion 
is  reduced.  Owing  to  the  lower  methane  conversion  at  higher 
Reynolds  numbers,  the  methane  CPO  tends  to  reduce  the  formation 
of  CO  and  H2,  as  observed  in  Fig.  14. 

4.  Conclusions 

The  chemical  kinetics  simultaneously  considering  methane 
combustion,  steam  reforming  and  carbon  dioxide  reforming  has 
been  successfully  modeled  to  account  for  the  reaction  character¬ 
istics  of  methane  catalytic  partial  oxidation  on  a  Rh  catalyst.  When 
the  Reynolds  number  is  fixed  (i.e.  Re  =  200),  the  methane  CPO  is 
activated  at  around  550  K  and  this  excited  temperature  is  hardly 
affected  by  the  C/O  ratio.  With  increasing  the  inlet  temperature  to 
the  activated  temperature  followed  by  decreasing  the  inlet  temper¬ 
ature,  it  is  worthy  of  note  that  different  paths  of  the  CH4  conversion 
and  the  maximum  temperature  are  exhibited,  yielding  a  hysteresis 
effect.  The  hysteresis  effect  developed  is  mainly  due  to  the  exother¬ 
mic  reaction  of  methane  combustion.  It  follows  that  the  historical 
trajectory  of  energy  state  in  the  catalyst  bed  plays  a  vital  role  in 
determining  the  result  of  methane  CPO.  When  the  inlet  temperature 
is  300 1<  and  C/O  =  1,  the  values  of  the  CFI4  conversion  in  the  upper 
and  lower  branches  are  64  and  0%,  respectively.  Therefore,  from  the 
viewpoint  of  methane  CPO  operation,  the  upper  branch  is  a  better 
condition  for  generating  syngas  and  hydrogen.  The  reaction  con¬ 
tours  indicate  that  methane  combustion  behaves  as  a  flame  sheet, 
which  takes  place  in  the  vicinity  of  the  reactor  entrance,  whereas 
steam  reforming  and  dry  reforming  prevail  over  the  entire  catalyst 
bed.  It  is  thus  recognized  that  the  methane  CPO  is  controlled  by 
the  combustion  reaction.  As  a  whole,  an  increase  in  inlet  temper¬ 
ature  facilitates  the  formation  of  FI2  and  CO.  This  arises  from  the 
fact  that  increasing  inlet  temperature  enhances  the  methane  com¬ 
bustion  and  thereby  increases  the  formation  of  H20  and  C02.  The 
increased  H20  and  C02  further  facilitate  the  generation  of  syngas, 
resulting  from  the  endothermic  reactions  of  steam  reforming  and 
carbon  dioxide  reforming.  Decreasing  C/O  ratio  renders  the  reac¬ 
tion  to  approach  complete  combustion  from  partial  oxidation;  as  a 
result,  the  hysteresis  effect  becomes  more  significant.  Alternatively, 
the  hysteresis  effect  withers  if  the  C/O  ratio  goes  up.  On  the  other 
hand,  increasing  the  Reynolds  number  results  in  the  decline  of  the 
hysteresis  effect  due  to  the  enhancement  of  the  inertial  force  of  the 
inlet  stream.  Flowever,  more  methane  is  depleted  and  more  heat  is 
released  as  the  Reynolds  number  increases,  the  hysteresis  behav¬ 
ior  of  maximum  temperature  is  thus  enlarged,  despite  the  lower 
methane  conversion. 
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